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ABSTRACT 


Relying  heavily  on  the  Elsasser  and  Culbertson  (1960)  computa- 
tional system,   computations  of  terrestrial  radiation  incident  at  a 
black  body  interface  have  been  programmed  for  an  arbitrary  atmosphe- 
ric sounding.   The  program  has  been  applied  to  a  random  sampling  of 
62  radiosondes  from  the  "Wark-sounding-catalog. " 

Since  Wark  et  al.  (1966)  have  kindly  made  available  simulated 
Nimbus  II  channel  2  and  channel  4  specific  intensities  for  these  same 
atmospheres,  a  multivariate  regression  was  derived  relating  downward 
flux  computations  to  Nimbus  II  channel  2  and  channel  4  readouts  as  well 
as  to  two  other  gross  air  parameters:   (a)  total  reduced  water  vapor 
depth,  and  (b)  interface  pressure.   All  four  independent  variables  gave 
high  statistical  significance,  with  channel  4  filtered  flux  accounting 
for  the  major  portion  of  the  "explained  variance"  of  the  dependent 
variable.   Total  ozone  was  also  tested  but  yielded  no  statistical  sig- 
nificance on  the  regression. 
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1.   Introduction 

A  knowledge  of  downward  infrared  radiation  flux  is  required  in 
order  to  understand  properly  the  heat  partitioning  in  the  earth-atmo- 
sphere and  sun  system.   The  atmosphere  is  much  more  transparent  to 
short-wave  than  to  long-wave  radiation.   This  is  depicted  by  the 
absorption  spectra  for  the  atmospheric  gases  on  p.  153  of  Fleagle  and 
Businger  [6].   From  these  spectra,  one  can  see  that  there  is  strong 
absorption  in  the  infrared  by  water  vapor,  carbon  dioxide  and  to  a 
lesser  degree  by  ozone.   It  is  the  goal  of  this  paper  to  develop  a 
computer  program  that  is  readily  adaptable  to  computation  of  downward 
radiation  flux  from  soundings  using  input  sounding  data  from  each  of 
these  three  atmospheric  constituents.   As  a  by-product,  the  secondary 
goal  is  to  relate  empirically  the  results  to  a  regression  equation  of 
selected  satellite  air  mass  properties. 

Wark  et  al.  [16]  studied  intensively  the  problem  of  upward  flux 
and  consequently  compiled  106  randomly  distributed  radiosonde  observa- 
tions from  the  I.G.Y.  files  in  what  is  hereby  termed  the  "Wark-catalog. " 
This  composite  set  of  soundings  cover  a  wide  range  of  radiative  condi- 
tions.  Moreover,  simulated  Nimbus  II  data  was  kindly  provided  by  J.H. 
Lienesch  of  the  National  Satellite  Center  for  the  same  106  soundings. 
Downward  infrared  flux  measurements  may  be  made  by  the  Kuhn-Suomi 
radiometer  but  these  give  only  limited  areal  coverage  compared  to  the 
MRR  data  from  satellites;  thus,  it  was  decided  to  compute  downward 
flux  in  a  manner  similar  (not  identical)  to  that  employed  by  Wark 
et  al.  in  his  1962  study,  and  to  perform  a  multiple  linear  regression 
analysis  with  certain  simulated  Nimbus  II  channel -readout  values  as 
input  variables.   The  computation  of  downward  infrared  flux  is  made 
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applicable  to  the  top  of  the  undercast  cloud  or  to  the  ground  level, 
but  the  computer  programming  is  left  flexible  so  as  to  apply  to  down- 
ward flux  at  an  arbitrary  number  of  sounding  levels  above  the  interface, 

2.   Nature  of  the  Data 

The  106  radiosonde  observations  in  the  Wark-Qatalog  are  given  as 
profile  listings  of  pressure  (mb) ,  temperature  (°K) ,  mixing  ratio  of 
water  vapor  (g/kg) ,  and  mixing  ratio  of  ozone  (cm  NTP  per  mb) .   Each 
sounding  or  case  atmosphere  is  identified  by  a  number,  the  name  of  the 
radiosonde  station,  and  the  date  and  time  of  the  sounding.   The  sky 
conditions  is  indicated  by  "clear"  or  by  the  height  (in  mb)  of  the 
top  of  the  undercast  cloud.   Temperatures  and  pressures  are  given  at 
approximately  25  +_  5  selected  levels  as  well  as  at  the  interface. 
The  interface  is  defined  for  this  study  as  being  the  earth's  surface 
for  clear  cases  or  the  top  surface  of  the  undercast  for  cloudy  cases. 
Total  (corrected)  water  vapor  and  ozone  masses  are  also  given  in  this 
Wark- catalog.   The  models  developed  by  the  Meterorological  Satellite 
Laboratory  were  taken  from  radiosonde  data  over  all  latitudes  and 
seasons.   All  the  soundings  were  from  the  I.G.Y  record- file  and  went 
to  at  least  25  mb.   Above  this  level,  temperatures  were  extrapolated 
parallel  to  the  appropriate  Supplementary  Standard  atmosphere  to  the 
pressure  of  0.1  mb.   Furthermore,  stratospheric  humidity  was  extrapo- 
lated to  conform  in  general  to  the  Supplementary  Standard  distribution 
for  the  particular  air  mass  at  the  station  (standard  cP,  standard  mT 
air  masses,  etc.).   The  62  cases  selected  are  composed  of  47  undercast 
cases  at  various  altitudes  and  15  clear  cases,  the  latter  being  ran- 
domly selected  from  the  Wark- catalog. 
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Nimbus  II  satellite  data  is  used  to  compute  filtered  channel  4 
upward  radiation  flux  and  the  interface  Black  Body  radiation  flux  for 
each  of  the  62  selected  atmospheres.   The  Nimbus  II  data  lists  total 
and  filtered  channel  4  (5y  to  30y)  infrared  specific  intensities  for 
five  zenith  angles  (theta) ,  that  is,  g  =  0° ,  20°,  45°,  60°,  and  78.5°. 
The  specific  intensities  for  9  =  90°  are  obtained  by  a  Lagrangian 
form  of  extrapolation  which  is  explained  below  in  connection  with 
equations  (19)  and  (20) .   The  radiation  intensities  of  the  62  cases  are 
listed  in  ergs/(cm2sec  stdn)  and  are  converted  to  watts/ (m2  stdn)  in 
this  paper.   It  should  be  noted  that  the  characteristics  of  the  sat- 
ellite infrared  sensing-filter  system  must  be  eliminated  so  that  the 
results  do  not  depend  upon  this  optical  system  to  any  appreciable  degree. 
The  National  Satellite  Laboratory  minimized  the  channel  2  and  channel 
4  sensor-filter  effects  on  the  Nimbus  II  data  by  normalizing  the  filter- 
sensor  response  functions  for  the  various  IR  wave  channels,  hence, 
the  resulting  response  was  the  "best  possible"  estimate  of  the  thermal 
flux  (5-30y)  for  the  channel  4  radiometer,  and  of  the  black  body  flux 
(10-lly)  or  its  equivalent  black  body  temperature  for  the  channel  2 
radiometer.   Also,  the  sounding  data  lists  temperature  at  the  inter- 
face of  each  of  the  case  atmospheres,  which  Keith  [9]  found  to  be  very 
accurately  described  by  the  Nimbus  II  channel  2  readout. 

3.   Computation  of  Downward  Flux  of  Radiation 

Downward  infrared  flux  of  radiation  arriving  at  the  interface  is 

quantitatively  computed  for  each  of  the  62  selected  atmospheres  from 

the  Wark-catalog  by  the  Elsasser  computational  system.   The  downward 

flux  computation  is  divided  into  five  parts . 
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Fd  =   F    +   Fco-   +   F03   -  AF(u*TVjt*)  -  AF(u*,tf)   (1) 
wv        * 

The  last  two  terms  of  equation  (1)  represent  the  "so  called  overlap 
of  flux"  and,  hence,  are  subtracted  from  the  first  three  terms  of 
equation  (1)   which  represents  the  downward  infrared  flux  not  cor- 
rected in  any  way  from  the  manner  described  by  Elsasser  and  Culbert- 
son  [5]  (henceforth  denoted  by  EC  for  abreviation) .   The  computation 
of  the  terms  of  equation  (1)  is   explained  below  in  connection  with 
equations  (6)  through  (11) . 

In  actuality,  an  energy-sharing  process  occurs  in  the  15u  and 
9.6p  bands  which  correspond  to  the  CO2-H2O  and  O3-H2O  overlap  fluxes, 
respectively,  thus,  accounting  for  the  AF-terms.   For  example, 
according  to  Hanel,  Bandeen,  and  Conrath  [7],  ozone  should  be  con- 
sidered the  primary  radiator  in  the  9.6y  band,  hence,  water  vapor 
is  the  weak  secondary  absorber-radiator  resulting  in  the  correction  due 
to  the  shared  radiation,  AF (u* , u*) • 

It  is  necessary  to  know  the  layer  values  of  reduced  optical  path 
for  water  vapor,  ozone,  and  carbon  dioxide  in  the  computation  of  radia- 
tive fluxes  by  each  of  the  three  constituents.   The  reduced  optical 
path  du*  of  an  absorber  has  been  taken  from  EC,  and  its  relation  to  the 
true  optical  path  is  given  by 

du*     =     P    faoVdu 

Po\T/  (2) 

According  to  EC  [5],  the  dependence  on  temperature  is  slight  and,  hence, 

is  excluded  in  this  study.   The  finite  difference  relations  for  the 

resulting  reduced  optical  paths  of  the  three  main  constituents  are  given 

by 

K       K 
u*  =  SAiU*i  =  E(qi+ii)(pi+i5)(Ai+i3P)  X  1.00705  X  10-6(g/cm2)     (3) 
i=l     i=l 
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K        K 
U*  =   ZAiU*  =   Z(Qi+jj)  (Ai+ilp)  X  10"5     (cm  NTP)  (A) 

i=l      i=l 


K        K 
V*   =   EA/U*  =   E(pi+JJ  (Ai+i5p)  X  2.4469  X  10-4   (cm  NTp)     (5) 
i=l      i=l    2 


Basically,  the  EC  computational  technique  for  the  determination  of 
downward  infrared  radiation  flux  is  realized  by  the  integration  of 
the  area  under  the  curve  as  shown  in  Fig.  1  with  the  use  of  Elsasser- 
Culbertson  [5]  flux  tables  of  R(u*,T)  which  have  been  incorporated 
in  the  author's  computer  program  RADFLUX  found  in  Appendix  I-A.   In 
essence,  one  enters  the  EC  tables  with  logarithmic  values  of  the 
constituent  optical  paths  and  their  corresponding  temperatures,  and 
obtains  the  respective  "R"  values  which  are  explained  below  in  connec- 
tion with  equations  (6)  through  (11)  .   These  "R"  values  are  then 
integrated  over  the  entire  temperature  range  of  the  sounding  result- 
ing in  the  flux  of  radiant  energy  due  to  the  particular  constituent. 
There  are  various  forms  that  can  be  used  to  determine  the  energy 
integral  but  Elsasser  et  al.  [5]  and  Yamamoto  [17]  used  temperature 
instead  of  optical  path  (u)  as  the  independent  variable. 
The  "R"  term  has  been  taken  from  EC  and  is  given  by 


R(u*,T)   -  f £  dB,      (u 


F(u*)d„  ,  (6) 

hence,  one  can  readily  see  from  equation  (6)  that  infrared  flux  can 
be  realized  by  an  integration  scheme  of  the  R(u*,T)  term  over  the 
temperature  range.   For  example,  downward  infrared  radiation  flux 
for  water  vapor  and  its  corresponding  finite  difference  form  are 
given  by 
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Fig.  1.   Schematic  form  of  the  distribution  of  R(U^T)  from  tables  in 
Elsasser  and  Culbertson  [5j  as  applied  to  a  typical  sounding  representative 
of  the  U-T  distribution  for  any  of  the  three  atmospheric  constituents. 
The  area  bounded  by  the  sounding  and  the  isopleths  U*  =  0  and  U*  =  Uf 
is  the  downward  flux  uncorrected  for  possible  overlap. 
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Fw  =  i  R(uf*,T)dT  +   3  °r£*(T),tJ  dT  (7) 

Tf 

i=if  i=80°C 

Fwv  ■  z  ^i-^^*±+^T±+h>   (Ti-Ti-i)  +  E  H+h^*f^±)^±+i-h   + 

i=l  i=if 

-2 73oC 

-8<JC  R(uf*,T)dT   (8) 

For  water  vapor,  the  range  of  wave  number  (\b/Vf)  is  listed,  by  40 

cm~l  intervals,  from  20  to  2600  cm"-'-,  which  is  regarded  as  spanning 

the  effective  infrared  range  of  the  terrestrial  radiation  spectrum. 

The  finite  difference  approximation  for  the  R-^  +  u  in  equation  (8) 

assumes  the  form  given  by 

Ri  +  Jj  =  [Ri(ui*,Ti)  +  Ri  +  1  (ui  +  i*,Ti  +  ]_)]  /2         (9) 
and  the  last  term  of  equation  (8)  is  obtained  from  Table  20  of  EC.   A 
similar  formulation  for  Fqqo  an<^  Foi  ■*-s  realized  by  the  use  of  Tables 
11  and  13,  respectively,  in  EC;  however,  Rco2  and  ^3  values  are  re- 
stricted to  the  frequency  ranges  of  540-820  cm-1  and  970-1130  cm   , 
respectively. 

The  last  two  terms  of  equation  (1) ,  the  H2O-CO2  and  H2O-O3  over- 
lap flux  terms,  are  defined  by  the  following  expressions 
Tf  _  TQ  _ 

'(11*01*)   =  J   TF(u*f)RC02(K*f,T)dT  +  J  xF(u*)RC02(U*,T)dT   (10) 

Tf 


AFI 


Tf=  T' 


AF(u*,U*)=  J   TF(u*f)R03(U*f,T)dT  +  J  t¥(u*)Rq3(V* ,T)dT  (11) 

Tf 

Moreover,  values  of  £  In*)    in  the  H2O-CO2  overlap  was  determined  by 
Martin  [10]  according  to 


Av  /       (12) 


TF(u*)  =  i  rFvBvTAv   =  tj^^f  "*"**   dt)  BvT 


vf 

BVTAv 


£   B,,t  Av 


v0   v 
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where  L   is  the  generalized  absorption  coefficient  listed  by  40  cm"^ 
intevals  and  B  m  is  the  Planck  specific  black  body  flux  for  a  given 
temperature.   The  results  of  this  computation  are  shown  in  Table  1. 
The  H2O-CO2  overlap  flux  occurs  in  the  540-820  cm   range,  hence  the 
water  vapor  TF(u)  within  this  range  was  obtained  by  performing  the 

summation  indicated  by  equation  (12)  using  both  T   ,  and  B   at  20° C, 

3       lFv       vT 

and  at  40  cm   intervals.   Lv  values  were  taken  from  Table  10  in  EC. 

According  to  EC,  20°C  is  a  representative  base  temperature  for  T 

Fv 

observations  since  the  H2O-CO2  overlap  is  in  practice  most  important 
in  the  lower  layers  of  the  atmosphere.   In  contrast,  Wark  et  al .  [16] 
used  Yamamoto's  [17]  tFv)0  values  which  were  based  upon  theoretical 
transition  probabilities;  however,  Palmer's  [14]  laboratory  measure- 
ments made  at  300°K  seem  to  give  good  agreement  with  Yamamoto's 
values  for  25  cnT^  intervals.   After  some  additional  smoothing,  it 
was  noted  that  the  quoted  L  values  over  40  cm-1  intervals  in  Table 
10  of  EC  corresponded  closely  to  Yamamoto's  quoted  values  of  "j,  /2" 
because  of  the  slightly  wider  interval  used  in  EC. 

Furthermore,  iptu*)  in  the  H2O-O3  overlap  is  given  by 
t~f(u*)   =   exp[-0.1167  u*]  (13) 

Equation  (13)  uses  the  continuum  (970-1130  cm-1)  absorption  coef- 
ficient of  water  vapor  proposed  by  Hanel  et  al.  [7].   In  both  the 
H2O-CO2  and  H2O-O3  overlap  regimes,  the  water-vapor  sharing  path  has 
been  increased  by  the  factor  5/3  in  order  to  transform  beam  trans- 
missivity  to  flux  transmissivity . 

The  results  of  the  foregoing  Elsasser  computational  scheme 
using  the  author's  computer  program  RADFLUX  on  the  62  selected 
case  atmospheres  from  the  Wark-catalog  are  shown  in  Table  2  along 
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Table  1.   Water  vapor  flux  transmissivities  in  the  H2O-CO2  overlap 
regime  (540-820  cm-1) . 


log  u*     xpCu*) 


log  u' 


TF(ll*) 


-6.0 

0.9999997150 

-2.7 

0.9994259 

-5.7 

0.9999994259 

-2.3 

0.998527 

-5.3 

0.999998527 

-2.0 

0.997150 

-5.0 

0.999997150 

-1.7 

0.994259 

-4.7 

0.999994259 

-1.3 

0.985665 

-4.3 

0.99998527 

-.10 

0.97128 

-4.0 

0.99997150 

-0.7 

0.94431 

-3.7 

0.99994259 

-0.3 

0.86680 

-3.3 

0.9998527 

0 

0.77328 

-3.0 

0.9997150 

+0.3 

0.67567 

+0.7 

0.54833 

21 


with  their  corresponding  case-number  identifiers  from  this  catalog. 
The  table  lists  uncorrected  fluxes  for  water  vapor,  ozone,  and  car- 
bon dioxide  together  with  the  CO2  and  ozone  overlaps .   The  CO2  and 
ozone  overlaps  are  equally  partitioned  arbitrarily  between  CO2  and 
water  vapor,  and  ozone  and  water  vapor,  respectively.   As  a  result, 
half  of  the  overlapped  flux  is  restored  to  CO2  and  ozone,  that  is, 
the  partitioned  overlap  fluxes  are  subtracted  from  their  respective 
uncorrected  constituent  fluxes  as  shown  below  in  connection  with 
equations  (14)  through  (16),  thus,  realizing  corrected  constituent 
fluxes . 

Fwv(CORR)   =   Fwv  -  AF(u*,U*)/2  -  AF(u*,V)/2  (14) 

FC02(C0RR)   =   FC02  _  AF(u*,M*)/2  (15) 

F03(CORR)   =   F03  -  AF(u*,U*)/2  (16) 

It  should  be  noted  that  all  the  fluxes  shown  in  Table  2  are  in  watts/ 
m  and  that  the  mean  of  each  overlap  and  constituent  flux  as  well 
as  the  average  percent  contribution  of  the  corrected  constituent 
flux  on  the  total  downward  flux  are  given  at  the  bottom  of  the  table . 
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H2U-C02 
OVLKLAP 


2 
3 
4 
7 
8 
10 
12 
13 
20 
23 
27 
31 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
ICO 
MEAN 
VALUE 


51 
47 
53 
33 
36 
57 
43 
^0 
54 
45 
34 
15 
4C 
37 
3? 
45 
32 
45 
24 
17 
32 
5S 
31 
43 
39 
47 
18 
^7 
?4 
30 
20 
26 
2  8 
29 
32 
40 
49 
55 
47 
31 
30 
24 
43 
23 
18 
36 


.423 
.r>97 
.  5  6  3 
.^■7f 
.3?  7 
.929 
.u63 
.274 
.531 
.347 
..A3 
.7^8 
.182 
.616 
.215 
.77^ 
.  Sr<4 
.9  3.'' 
.351 
.377 
.795 
.728 
.  L: .'  5 
.924 
.92  7 
.  3 44 
.276 
.56C 
.158 
.741 
.222 
.922 
.2^9 
.66  3 
.466 
."31 
.5  66 
.726 
.391 
.5  15 
.767 
.522 
.949 
.776 
.2^3 
.487 


36.63" 
35.667 
24.369 


12. 
24 

24 
23 

63 
45 


39C 
563 
829 
987 
733 
5Z1 


47.421 
53.  147 

.265 
488 


42 

50 
41 
43 
42 
32 


261 
122 
395 
643 


H20-U3 
OVERLAP 

3.066 
4.138 
4.6C1 
1.395 
3.652 
5.143 
3.861 
5.956 
7.672 
5.519 
4.074 
1.294 
3.P46 
3.922 
4.797 
6.5  30 
2.2  26 
5.642 
1.0  99 

1.5  64 

3.6  40 
7.619 
4.055 
5.251 
4.713 
6.044 
1.816 
5.:  57 
3.213 
2.fc78 
1.615 
2.752 
3.C66 
1.939 
4.020 
3.324 
5.531 
5.573 
5.332 
5.311 
5.159 
4.796 
7.7C4 
2.H22 
1.166 

4.  146 
4.4C4 
3.:  66 

.261 

.963 

.177 
2*855 
3.146 
5.464 
5.780 
6 .  J  5  3 
7.433 
5.15J 
7.567 
1.2  59 
6  .  v  1  • ' 

5.  143 
5.   14  3 


UNCOR- 

9LCTE0 
HV  FLUX 
255. P<>8 
246.671 
168.837 
155.552 


123 
263 
l^C 
IriB. 


25  6 
871 
5-4 
8  31 


169.774 
136.574 
131.287 

71.  5S*. 
125.149 


4S, 

76 

12P. 

1  6 

2  7 

6.3, 

6P 
t.7 

214. 


4(>6 
86t 
7^9 
971 
4^7 
715 
e77 
821 
74  0 


78.76t 
145.731 


'J  ' 

1^6 
116 

Ht> 

44 

6  0 

M 

94 

66 
142 

2   r. 

197 

lb6.667 
7i  .492 
f.6.361 
'1  3  .  1 9  1 

121. 
48 
3  5 

llr 
92 

1     1 


916 
71b 
113 
2^t 
451 
615 

on 

*_   ~ 

17C 

016 
833 
163 
664 
2'3  7 


111 
124 
946 
217 
464 
135 
59.652 
17.517 


6  7 

r-3, 
49, 

2  01, 

147 

14 

197, 


3  24 
94  9 
17? 
636 
232 
Cbb 
972 


136.471 
132.36,2 
177.674 
1 2  8  .  2  "•  3 
139.56  3 
66  .844 


CORRfCT- 
tn  WV 

FLUX 
228.653 
2  20.8~<4 
139.655 
138.  31  6 

1  ,2.26  9 

2  3  2.335 
1  6  7  .  34  2 
130.716 
158.672 
111.142 
1  12.229 

6  3.  >5  3 

1.3.135 
77.627 
58.362 

1'.  2.644 
69.43C 

161.671 
53.990 
59.4"7 
4  9  .  6  '  3 

182.966 
6  ".98 7 

171.14  3 


j.  -  J 

124 

20 

84 

99 

7: 


,598 

,J24 

,^67 

95  7 

765 

j5 


37.423 


4.1L7 


33.172 

45.392 

45.492 

78.1" 5 

6  8.58  8 

123.  5. '5 

173.126 

166.637 

1  6  3  .  5 •:  6 

52.    76 

48.39e 

25.532 

95.28  5 

34.625 

26.215 

94. 9"- C 

71.^47 

81.759 

47.33  7 

10.631 

54.654 

4..  117 

35.676 

226.239 

121.581 

113.329 

167.682 

1  14.754 

1 33.335 

156. 4"3 

1  "3.63  7 

1  15.794 

4  8.^0  1 

1^0.429 


121.194 
AVG    VALUc    CONTRIBUTION    i)  i    TOTAL 

DOWNWARD    FLUX     IN    PFRCfciT 73.97 


UNCOR- 
RECTED 

oFflux 

10.874 

11.410 

9.537 

8.094 

7.774 

13.398 

11.503 

10.081 

11.994 

9.247 

9.325 

6.281 

?:!H 

7.832 

9.811 

7.264 

13.859 

5.794 

7.603 

5.630 

13.182 

7.678 

10.663 

10.050 

11.364 

5.194 

8.305 

9.281 

8.055 

5.850 

6.394 

6.678 

7.559 

7.643 

8.855 

10.767 

11.145 

10.664 

7.586 

6.794 

5.6  04 

9.886 

6.  196 

5.228 

8.9C4 

8.C34 

7.638 

5.535 

3.274 

5.931 

6.392 

6.354 

11.412 

10.253 

9.121 

12.978 

10.134 

13.082 

10.920 

9.542 

10.008 

6.553 

8.779 


CORRECT' 
ED  03 
FLUX 
9.341 
9.341 
7.136 
7.397 
5.948 

1'  .826 
9.573 
7.103 
8.158 
6.487 
7.289 
5.634 
7.908 
7.31C 
5.434 
6.546 
6.150 

11.038 
5.244 
6.821 
3.8i: 
9.273 
5.651 
8.037 
7.694 
8.342 
4.287 
5.777 
7.675 
6.617 
5.042 
5.018 
5.145 
6.590 
5.633 
7.193 
P.C01 
6.359 
7.99C 
4.931 
4.214 
3.206 
6.034 
4.785 
4.644 
6.830 
5.832 
6.095 
5.405 
2.782 
5.842 
4.964 
4.781 
8.680 
7.363 
6.094 
9.261 
7.559 
9.298 

13.290 
6.537 
7.436 
4.C31 

6.726 

5.96 


UNCOR- 
RECTED 
CC2  FLUX 
73.431 
71.421 
56.286 
56.826 
48.5  50 
69.619 
62.955 
54.728 
61.467 
51.151 
52.046 
32.212 
48.695 
45.701 
31.682 
51.199 
42.632 
65.681 
30.912 
27.909 
26.384 
61.622 
31.367 
48.377 
46.979 
56.187 
19.450 
37.569 
43.899 
35.407 
22.546 
26.2C  6 
27.723 
36.166 
35.273 
46.774 
61.434 
63.o74 
61.123 
29.848 
27.835 
18.050 
46.-74 
24.261 
18.494 
46.103 
37.709 
38.045 
26.687 
14.466 
32.756 
25.524 
24.777 
73.731 
48.342 
46.967 
61.539 
46.911 
59.655 
57.266 
48.626 
54.564 
26.749 

43.780 


C?8R!SJ- 

FLUX 

47.719 

47.622 

29.504 

4'  .287 

29.386 

40.6  55 

41.423 

29.591 

34.202 

28.478 

35.025 

24.356 

28.604 

26.793 

15.574 

28.31. 

26.205 

42.716 

18.737 

19.22; 

9.686 
33.756 
15.615 
26.415 
27.016 
32.515 
10.312 
18.789 
26.82 
20.037 
12.435 
12.745 
13.579 
21.225 
19.039 
26.756 
36.641 
35.81") 
37.423 
14.09 
12.45i 

5.78^ 
24.C99 
12.374 

9.347 
27.859 
19.394 
20.211 
14.502 

8.271 
23.477 
13.11  i 
12.783 
41.864 
25.582 
23.256 
34.966 
25.768 
34.411 
36.626 
27.265 
33.366 
10.428 

25.069 

23.07 


TOTAL 

DONNWARO 
FLUX 
285.713 
277.767 
176.296 
185.999 
137. oC2 
283.815 
218.338 
167.410 
2  01.332 
146.107 
154.542 
93.  -46 
139.647 

1  1 1 .  7  30 
79.370 

137.500 

121.785 

235.425 

74.971 

85.448 

63.3  99 

225.997 

82.252 

155.596 

135. 3C7 

164. a81 

34.666 

109.523 

134.260 

96.659 

50.650 

63.154 

64.217 

105.920 

93.260 

154.456 

2  17.768 
210.805 
205.932 

71.100 

65.  63 

37.526 

125.419 

51.983 

40.207 

129.589 

97.174 

108.365 

67.244 

21.8  84 

80.973 

5  8.180 

53.170 

276.784 

154.526 

142.679 

211.909 

148. .81 

147.  45 

203.319 

137.440 

156.596 

62.46C 

132.223 


Table  2.    Selected  62  caaa  atmospheres  as  indantiflad  from  the  Wark-catalog ,  and  the  computed 
downward  flux  from  program  RAD  FLUX. 
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4.   Determination  of  Channel  4  Upward  Flux  of  Radiation. 

For  each  sounding  case,  the  filter-normalized  upward  flux  in  the 
wavelength  of  5y  to  30u  was  available  as  a  result  of  computations  fur- 
nished by  Lienesch.   This  was  in  the  form  of  the  "so-called"  channel 
4  specific  intensity  at  the  top  of  the  atmosphere.   Channel  4  specific 
intensity  14(6)  was  then  transformed  to  a  channel  4  flux  for  each  of 
the  selected  62  cases  by  using  the  trapezoidal  rule  for  integration  on 

14(8)  as  given  by 

-  2 

4>F4   =   2tt      /2   I4(6)    d(sin   6)  (17) 

o  2 

and  its  corresponding  finite  difference  form  is  given  by 

<j>f4  =  2il  i    i4(e.)+i4(ei+i) 

|ei=0 - Xis(sinz{ei+1}    -   sin^Oi})]  (18) 

where  I  =  1,  2,  3,^.^5,  and  6i  =  0°,  62  =  20°,  63  =  45°,  94  =  60°, 
85  =  78.5°,  and  6^  =  90°.   The  trapezoidal  rule  was  selected  instead 
of  other  finite  integration  schemes  because  of  the  unequal  zenith- 
angle  intervals  of  integration.   Also  the  channel  4  upward  flux  relation 
is  the  same  as  that  used  by  Wark  et  al.   L.1§J  which  allows  for  the  usual 
approximation  of  azimuthal  symmetry. 

The  zenith  angles  (6)  of  the  Nimbus  II  listings  went  as  far  as 
78.5°,  hence,  the  specific  intensity  f or  0  =  90°  was  extrapolated  by  a 
Lagrangian  form  expressed  by 

i(e)  =  ii(e-e2)(e-e3)  +  i2(e-e1)(e-e3)  +  13(6-01) (e-e2)      (19) 


(61-62)01-63)     (62-61)  (62-63)    (63-61)  (63-6-j) 
and  which  reduces  to 

1(90°)  =  I78.50  +  0.6866l45o  -  1.8649160°  +  1-1783l78.5°        (20) 
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According  to  Wark  et  al.  ^16j ,  radiation  for  angles  greater  than 
78.5°  contributes  only  about  4  percent  of  the  total  flux  and  the  energy 
obtained  by  integrating  eq.  (17)  between  0  =  78.5°  and  6  =  90°  lies 
well  within  this  4  percent  range. 

The  black  body  radiation  flux  (oTg)  at  the  interface  is  given 
directly  by  the  interface  temperature  (°K) ,  and  when  o=   0.56687  X  10"^ 
(watts/m2  °K  )  ,  oTg  is  in  watts/m2. 

The  two  satellite  air  mass  properties,  <f>F,  and  aT™,  are  made  avail- 
able for  this  study  by  the  foregoing  procedures,  where  Tg  is  taken 
equivalent  to  the  channel  2  emittance  temperature  in  a  simulated 
Nimbus  readout.   Moreover,  the  foregoing  integration  employing  eqs. 
(18)  and  (20)  were  incorporated  in  the  author's  computer  program  SUMFLUX 
(Appendix  I  C) .   The  62  selected  Nimbus  II  cases  as  well  as  the  computed 
air  mass  properties  are  shown  in  Table  3. 

5.   Statistical  Results  and  Inferences 

The  object  in  the  statistical  treatment  which  follows  is  to  obtain 
a  specification  equation  for  F^  in  terms  of  satellite  readouts  and 
other  readily  accessible  air-mass  data.   To  accomplish  this,  there  was 
available  a  BMD03R  computer  program  that  executes  the  multiple  regres- 
sion analysis  establishing  a  regression  plane  given  by  the  linear  form 
of 

Y  =  AQ  +  A-^  +  A2Z2  +  A3Z3  +  A4Z4  +  A5Z5  (21) 

with  the  final  results  of  the  regression  plane  shown  in  Table  4. 

The  original  set  of  62  independent  data  atmospheres  was  split 
randomly  into  two  equal  sets.   One  set,  Sample  A  which  is  denoted  by 
the  asterisks  in  Table  3,  was  used  to  develop  the  regression  while  the 
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Table  3.      Channel  4  filtered  flux  and  interface  black  body  flux  on 
the  62  case  simulated  Nimbus  II  data. 

CASE  INTERFACE          INTERFACE           CHAN^FL    4  RLACK 

PRESSURE          TEMPERATURE  FLUX  BODY 

FLUX 

*  2  1CC0  289  181.076  395.435 

*  3  1009  289  186.376  395.435 

*  4  1014  284  171.133  368.771 

*  7  998  267  153.014  288.090 

*  8  1000  262  143.171  267.1C9 

*  10  1000  298  189.664  447.042 

*  12  923  282  170.847  358.492 

*  13  1003  278  163.286  338.581 
20  850  291  188.125  406.495 
23  1000  274  155.662  319.511 
27  941  258  140.013  251.167 
31  1020  235  105.440  172.884 

50  850  275  164.415  324.201 

51  850  268  156.219  292.430 

52  500  261  144.776  263.054 

53  830  275  163.762  324.2^1 

54  703  264  14C.085  275.359 

55  908  285  184.727  373.993 

56  526  249  120.077  217.912 

57  10C6  226  99.153  147.883 

58  400  257  133.449  2*7.295 

59  850  291  181.957  406.495 

60  500  261  141.883  263.054 

61  7C0  283  175.614  363.604 

62  652  274  162.747  319.511 

63  754  283  178.669  363.604 

64  400  230  96.914  1^8.633 

65  7C0  267  143.358  288.090 

66  700  266  145.276  283.798 

67  568  261  140.226  263.054 

68  500  23«  106.511  181.883 

69  7C0  248  119.152  214.432 

70  700  250  122.381  221.434 

71  476  264  146.354  275.359 

72  466  265  151.586  279.555 

73  7C0  280  165.731  348.430 

74  810  285  176.863  373.993 
■75  930  289  187.088  395.435 
76  932  276  164.492  328.942 

*  77  5C0  256  138.983  243.469 

*  78  5CO  253  131.614  232.255 

*  79  4C0  244  118.779  200.929 

*  80  8C6  268  |?f-??l  ?2£'5?§ 
•81  5C0  242  114.152  194.422 

*  82  400  233  101.118  167.073 

*83  800  264  \^A'\ll  ?tH?2 

*  84  722  263  143. '^31  271.210 

*  85  7C0  266  145.922  283.798 
*86  720  248  116.863  214.432 

*  R7  AGO  221  79.778  135.224 

*  8H  818  248  117.937  214.432 

*  §9  500  246  121.593  207.598 

*  90  370  251  130.516  224.998 
•92  7C0  284  178.260  368.771 

*  11  7C0  283  175.544  363.604 
•94  850  289  183.440  395.435 

*  95  700  278  167. 0C3  338.581 
*"  850  285  180.0C8  373.993 

m  lis  I!!  :  '  £*'.& 

*99  900  271  158.044  1^1$$ 

•100  700  255  132.669  239.687 
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remaining  set,  Sample  B,  was  used  to  test  the  stability  of  the  test 
data  on  the  check-data  which  was  also  derived  by  the  author's  computer 
program  RAD FLUX. 

According  to  Miller  \_13j  and  Martin  et  al.  JJ.1J  ,  the  regression 
screening  analysis  uses  a  F1  critical  value  both  to  accept  and  remove 
individual  variables  from  the  regression.   The  critical  F1  value  uti- 
lizes a  "probability  level"  that  is  a  function  of  the  number  "k"  of 
predictors  tested  and  of  the  Sth  selected  prediction  system.   The  pro- 
posed critical  value  corresponds  to  the  tabulated  F'  value  at  a  prob- 
ability level  of   a=a*/(k-S+l)  with  1  and  N-S-l  degrees  of  freedom, 
where  N  is  the  sample  size  to  be  tested.  Since  the  predictors  (k=5) 
have  been  pre-selected,  the  selection  choice  "S"  is  equal  to  1,  hence 
K-S+l  is  equal  to  5. 

The  predictors  employed  are  <J>F4,  oTg,  u*,  the  interface  pressure 
"p",  and  U*  which  correspond  to  Z]_,  Z£,  Z3,  Z4  and  Z5,  respectively,  in 
eq.  (21).   The  order  of  the  predictors  is  a  vital  consideration  in 
describing  the  "explained  sum  of  squares  by  the  regression."   Channel 
4  filtered  flux  (tyF^)    is  referred  to  as  the  thermal  channel,  for  its 
infrared  range  is  5u-  30y,  whereas  channel  2  filtered  flux,  the  window 
channel,  has  an  infrared  range  of  lO-llp.   Channel  2  filtered  flux, 
however,  is  equivalent  to  the  interface  black  body  flux,  aTg  ,and  there- 
fore contains  very  little  atmospheric  structural  information.   Pre- 
scribing <(>F4  ahead  of  oTg  allows  (JjF^  to  assume  a  larger  proportion  of 
the  explained  variance  than  aTg. 

The  most  readily  adaptable  test  for  the  significance  of  the 
predictor  is  the  F  test  based  upon  an  analysis  of  variance.   The  F 
statistic  is  defined  as 
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F  =  Mean  squares  explained  by  the  predictor (s) (22) 

Mean  squares  not  explained  by  the  predictor (s) 

The  F  statistics  derived  from  the  results  of  the  BMD03R  program  are 

shown  in  Tables  4  and  5.   These  statistics  are  used  to  determine  the 

significant  specification  on  F^  by  the  multivariate  regression  and  by 

each  of  the  predictors.   Significance  in  this  analysis  means  that  the 

independent  variable  and/or  variables  selected  for  the  F  test  did  not 

affect  the  outcome  of  the  regression  by  mere  chance  at  the  specified 

level  of  belief.   For  this  analysis,  an  "a"  of  0.01  is  chosen  resulting 

in  a  composite  99%  confidence  level  of  belief  for  the  multiple  regression, 

Also  shown  in  Tables  4  and  5  are  the  tabulated  critical  F'  values 

and  the  "a"  critical  level  for  rejection  of  a  possible  air-mass  variable 

together  with  the  corresponding  "F-value  upon  entry."   Since  "a*" 

is  chosen  to  be  0.01,  the  "a"  critical  levels  for  the  5  predictor 

regression  are  0.002.   The  "F-values  upon  entry"  are  defined  by 

F(l,N-m-l)   =  Mean  Squares  explained  by  the  regression (23) 

Mean  Squares  not  explained  by  the  regression 

where  "N"  denotes  the  number  of  cases  in  the  sample,  "m"  is  the 

sequential  order  of  the  predictors  (m  =  1,  2....k),  "k"  is  the  number 

of  predictors  in  the  regression,  and  the  number  pair  (l,N-m-l)  denotes 

the  number  of  degrees  of  freedom  associated  with  the  numerator  and  the 

denominator,  respectively.   For  example,  the  "F-value  upon  entry"  for 

the  first  predictor  (<}>Fa)  in  the  regression  is  given  by 

(24) 
F(l,29)   =   (Sum  of  Squares  attributable  to  the  regression)/! 

(Sum  of  Squares  not  attributable  to  the  regression) /29 

The  addition  of  the  next  predictor ,oTg,  requires  the  definition  of 

4 
a  new  F-statistic  specifically  restricted  to  the  testing  of  aTg, 

after  the  contribution  of  (JjF^  has  been  excluded.   This  new  F-statistic 

follows  the  same  form  as  given  by  eqs.  (22)  and  (23)  but  with  one  less 
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degree  of  freedom  in  the  denominator.   The  remaining  independent 

variables  are  examined  successively  in  the  same  manner  with  a  new 

F-statistic  specifically  restricted  to  the  testing  of  the  selected 

variable  after  the  contribution  of  the  variables  already  tested  have 

been  excluded. 

The  composite  F-statistic  for  the  multiple  regression  is  based 

upon 

F(k,N-k-l)  =  (Sum  of  Squares  attributable  to  the  regression) /k ' 

(Sum  of  Squares  not  attributable  to  the  regression) / (N-k-1) 

Use  of  eq.  (25)  with  the  results  from  Table  4  gives  the  F-statistic 

for  the  5  predictor  multiple  regression  as  F  =  116.03  which  is  far  in 

excess  of  a  critical  F'(5,25)  =  3.86  appropriate  at  a  99%  confidence 

level. 

The  results  of  Table  4  show  not  only  that  the  multivariate 
regression,  but  also  that  each  of  the  single  variate  predictors  are 
significant  with  the  exception  of  the  last  independent  variable,  U*. 
The  "F-statistic  upon  entry"  for  u*  is  0.90416  in  contrast  to  the 
required  (tabulated)  critical  F' -value  of  14.0  at  99.8%  confidence  level 
of  belief.   This  indicates  that  U*  could  only  have  affected  the  outcome 
of  the  multivariate  regression  by  chance  at  a  very  low  level  of  probab- 
ility.  An  explanation  for  the  lack  of  significance  by  ozone  could  lie 
in  that  the  computed  amounts  of  ozone  in  a  column  of  atmosphere  are  near- 
ly constant  across  the  62  case  sample  resulting  in  a  negligible  influence 
on  the  regression.   Since  U*  showed  so  little  significance  in  the 
regression  analysis,  another  analysis  was  made  with  U*  excluded.   The 
regression  plane  for  the  4  predictor  case  is  given  by  the  linear 
expression 

Y'   =  A0  +  AjZi  +  A^Z2  +  A^Z3  +  A^Z4  (26) 
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with  the  final  results  of  the  multivariate  regression  shown  in 
Table  5.   The  same  procedures  and  testing  criteria  as  before  were 
used  resulting  in  each  of  the  4  predictors  as  well  as  the  composite 
set  of  predictors  showing  significant  contribution  on  the  regression, 
thus,  the  4  predictor  regression  was  selected  as  the  testing  regres- 
sion plane  for  the  independent  data  sample.   The  five  variate  regres- 
sion equation  is  given  by 

Fd  =  -85.535  +  1.750(<}>F4)  -  0.486(ct£)  +  47.466(u*)  +  0.129(p)  + 

-121.866(U*)  (27) 

and  the  four  variate  regression  equation  is  given  by 

Fd  =  -128.089  +  2.022((J>F4)  -  0.605(ot£)  +  49.897(u*)  +  0.129(p)  (28) 
It  is  interesting  to  note  that  the  multiple  correlation  coefficients 
shown  in  Tables  4  and  5  decrease  very  slightly  from  0.9791  in  the  five 
variate  case  to  0.9784  in  the  four  variate  case.   An  explanation  for 
this  could  lie  according  to  Panofsky  et  al.  [l5J  and  Martin  et  al.  [llj 
that  the  longest  equation  may  have  actually  overfitted  the  first  analysis 
ascribing  some  of  the  variation  due  to  small-scale  perturbation  to  U* 
by  chance. 

Finally,  the  4  predictor  regression  was  tested  using  Sample  B 
data  by  an  analysis  of  variance  or  of  mean  squares,  that  is,  a  test 
of  variation  about  the  regression  plane.   Again,  significance  of  the 
regression  plane  was  tested  by  the  F  test  as  given  by  equation  (24) . 

Total  sum  of  squares  with  its  components  is  given  by 

N      _       N  N      _ 

E  (y,-y)2  =  E  (y.-y'.)2  +  E  (y.-y)2  (29) 

1-1  x        i=l  x  x      1=1  x 

(A)  (B)  (C) 

where  "y'."  denotes  the  regression  estimates  for  Fj,"y^"  denotes  the 
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observed  F^  values  as  determined  by  the  computer  program  RADFLUX,  and  "N" 
denotes  the  total  number  of  cases  in  the  independent  data  sample.   Term 
(A)  gives  the  total  sum  of  squares;  term  (B)  gives  the  sum  of  squares 
of  deviations  from  the  regression  plane  or  the  sum  of  squares  not 
explained  by  the  regression;  and,  term  (C)  gives  the  sum  of  squares  of 
deviations  of  the  regression  estimates  from  the  mean  (y)  of  the  observed 
values  or  the  sum  of  squares  explained  by  the  regression. 

The  statistics  on  the  independent  data  sample  are  given  in  Table 
6.   The  required  (tabulated)  critical  F'value  LF'(4,262/is  equal  to  4.14 
at  99%  confidence  level  of  belief.   One  can  readily  see  that  the  "F" 
statistic  for  the  independent  data  F(4,26)  =  118.2546,  is  considerably 
larger  than  its  corresponding  tabulated  critical  F'value,  thus  the  4 
predictor  regression  verifies  as  being  significant  for  determining 
Fj  within  99%  confidence  level  of  belief. 

Furthermore,  the  results  of  the  4  predictor  regression  on  the 

independent  data  sample  were  tested  for  shrinkage  using  the  relation 

given  by 

R^  =  1  -  (Sum  of  Squares  not  explained  by  the  regression)     (30) 

(Total  Sum  of  Squares) 

"  2 
where  R^   is  the  effective  coefficient  of  determination  and  R^  is  termed 

the  "effective  multiple  correlation  coefficient."   The  results  are  shown 

in  Table  6. 

The  percent  shrinkage  between  the  explained  sum  of  squares  of 

2*2  2 

samples  A  and  B  is  given  by  R^  -  R^  .   Here  R/   is  the  coefficient  of 

determination  for  eq.  (28)  (Sample  A  data).   Table  6  shows  that  the 
regression  performed  exceptionally  well  on  the  independent  data  sample 
with  only  a  5.62%  amount  of  variance  shrinkage  and  with  a  R4  and  R4 
of  90.1%  and  95.1%,  respectively. 
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Table  6.   Analysis  of  variance  on  the  independent  data  (Sample  B) 


Source  of  Variation 


Degrees    Sum 

of       of 
Freedom Squares 


Variance 


Explained  Variance 


86065.300    2868.8430 


Unexplained  Variance    26 


9462.072     157.6892 


Total  Variance 


30 


95527.372    3026.5322 


F-value 
F(4,26) 


Critical 
F'  value 
F'(4,26) 


Coeff.  of 
Determina- 
tion (R42) 


Multiple      Percent 

Correl. 

Coeff.  (R4)    Shrinkage 


118.2546   4.14 


0.901 


0.951 


5.62 
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The  foregoing  analysis  showed  that  the  4  predictor  regression  was 
indeed  verified  on  the  independent  data  sample;  but,  the  process 
of  randomly  dividing  the  original  62  Nimbus  II  cases  in  order  to  realize 
this  verification  excluded  31  perfectly  good  cases  from  the  regression 
analysis,  thus,  it  was  decided  to  obtain  a  regression  based  on  the 
entire  62  Nimbus  II  cases.   One  should  anticipate  a  more  representative 
regression  than  the  one  previously  derived,  because  the  sample  size  is 
twice  as  large  as  that  used  in  the  original  analysis.   Although 
the  new  analysis  used  the  same  methods  and  testing  criteria  as  that 
used  in  the  31  case  sample,  it  was,  however,  confined  to  the  4  predictor 
analysis  because  of  the  earlier  findings  in  the  study.   The  final  regres- 
sion is  given  by 

F^  =  -129.847  +  1.824(<}>F4)  -  0.391(oTg)  +  43.293(u*)  +  0.097(p)  (31) 
and  its  corresponding  statistics  are  given  in  Table  7.   It  is  interest- 
ing to  note  that  the  multiple  correlation  coefficient  decreased  slightly 
in  going  from  the  31  case  analysis  to  the  62  case  analysis.   An  explan- 
ation for  this  must  lie  in  that  the  sample  B  "explained  variance"  was 
slightly  smaller  than  that  for  sample  A. 

6.   Conclusions 

The  programmed  computational  system  for  downward  infrared  radiation 
flux  due  to  the  absorbing/emitting  properties  of  water  vapor,  ozone  and 
carbon  dioxide  is  given  in  Appendix  I  A  by  the  author  s  computer  program 
RADFLUX.   This  program  may  be  readily  used  for  the  determination  of  Fj . 
The  only  necessary  inputs  are  pressure  (mb) ,  temperature  (°K) ,  water 
vapor  mixing  ratio  (g/kg) ,  and  ozone  mixing  ratio  (cm  NTP  per  mb)  for 
selected  levels  in  a  sounding.   Moreover,  this  program  can' be  easily 
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modified  as  shown  in  Appendix  I  B  in  order  to  determine  F  from  any  pre- 
selected level  to  the  top  of  the  atmosphere  (0.1  mb) .   It  should  be 
noted  that  the  mean  ratio  of  downward  radiation  flux  due  to  water  vapor 

(F  )  with  and  without  overlap  corrections  to  the  total  downward 
wv 

radiation  flux  (F.)  was  73.97%  and  85.6%,  respectively. 

Furthermore,  total  and  filtered  channel  4  upward  radiation  flux 
can  be  determined  by  using  the  author's  computer  program  SUMFLUX  found 
in  Appendix  I  C.   Also,  the  interface  black  body  flux  or  its  channel  2 
equivalent  temperature  can  be  readily  determined  using  this  same  pro- 
gram. 

Finally,  the  multivariate  regression  equation  developed  using 
Sample  A  Nimbus  II  air-mass  properties  for  the  determination  of  F, 
had  an  accuracy  of  95.1%  on  the  independent  test  data,  Sample  B.   This 
can  be  seen  in  the  foregoing  analyses.   The  explained  standard  devia- 
tion of  the  Sample  A  regression  was  95.84%,  but  the  resulting  explained 
standard  deviation  on  the  independent  sample  was  95.1%,  giving  a  loss 
of  accuracy  of  2.74%.   A  partial  explanation  for  the  loss  of  accuracy 
lies  in  the  process  itself  of  developing  a  regression  on  a  finite  ran- 
dom sample.   One  should  not  anticipate  a  regression  determined  from  a 
finite  sample  to  perform  exactly  the  same  way  on  another  sample  of  the 
same  population.   This  is  inherent  in  the  randomness  of  the  selection 
process  and  of  the  variations  within  the  population.   This  effect  can 
be  anticipated  to  be  partly  minimized  by  using  the  largest  possible 
sample.   Another  partial  explanation  could  lie  in  the  assumption  that 
the  regression  is  linear,  but  in  actuality,  F,  might  vary  non-linear ly 
with  one  or  any  combination  of  its  predictors.   This  would  indicate 
that  the  true  regression  is  some  sort  of  surface  and  not  a  plane, 
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however,  since  the  shrinkage  of  the  multiple  correlation  was  only 
2.74%,  the  regression  plane  appears  to  be  a  good  estimate  of  the  true 
surface  or  of  physical  reality.   With  the  advent  of  improved  satellite 
sensing  systems,  the  regression  equations  can  be  expected  to  be 
improved  by  an  analysis  similar  to  the  one  given  in  this  paper.. 
Furthermore,  an  improvement  on  the  regression  plane  may  be  anticipated, 
if  the  satellite  data  as  well  as  the  radiosonde  data  could  be  strati- 
fied to  a  common  level,  namely  that  of  the  interface-   This,  of  course, 
would  exclude  pressure  as  a  possible  air-mass  predictor,  however,  one 
could  reasonably  expect  this  stratification  to  reduce  the  scatter  about 
the  regression  plane.   Moreover,  since  the  effect  of  ozone  in  a  column 
of  atmosphere  gives  a  very  nearly  constant  contribution  to  F,  and, 
thus,  has  very  little  effect  on  the  regression  plane,  and  since  one 
should  expect  that  the  regression  based  on  the  largest  possible  sample 
to  be  the  most  representative  to  the  true  regression  surface,  the 
best  regression  for  the  determination  of  F,  is,  therefore,  the  4 
predictor  regression  based  upon  the  62  case  Nimbus  sample,  eq,.(31). 
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APPENDIX  I 

Significant  computer  programs  developed  for  this  study  (Sections 
A,  B,  and  C) . 
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Appendix  I,  Section  A.  Computer  Program  RADFLUX  (FORTRAN  63) 


-COOP. BOX  T,  TUPAZ  J .S/ 1S/2S .60 .40000. 
-FTN.L.E. 

PROGRAM  RADFLUX 
C      PURPOSE  OF  THIS  STUDY  IS  TO  EXAMINE  DOWNWARD  FLUX 

OD I  MENS  I  ON  TT(13)  .U(21)  .UU(16) ,UUU(23) .WATER ( 2 1 » 13 ) . OZONE ( 16.13) ♦ 
1C02(23.13) .RWATERI21) .R0Z116) ,RC02(23) »P(40) .PMEAN( 40 ) .PD I FF (40 ) . 
2T(4  ) ,Q(40) .QMEAN(40) » QDI FF ( 40 ) »CQ ( 40 ) .CQMEAN(40) »CODIFF(40) . 
3SFW(40) .SFOI40) .SFC(40 ) .CFW(40) .CFO(40) .CFC(40) ♦ FUT ( 40 ) » GUT ( 40 ) . 
4DUT(40) »FF( 1) »GG( 1 ) .DD(1) .COTAUF(21) .03TAUF(40) .RCOLAP( 1). 
5CO2LAP(21).O3RLAP(40) ,R03LAP(1 ) .UC(21) »COXMIS(21) .TCOLAP( 1 ) . 
6T03LAP( 1) 

DIMENSION  <ASE(40) .LEVEL(40) 

DIMENSION  Z<40) .S<40) .V(40)»R(40) 
OD I  MENS  I  ON  RH20(40) .RO3(40) .RRCO < 40 ) .TEMPI ( 40 ) . TEMP2 ( 40 ) .TEMP3(40) , 
1  ROPE (13) »SOPE( 13).TOPE(13).TTW(13).TTO(13).TTC(13).RW(13).RO(13). 
2RC(13) »RURF( 13) .SURF ( 1 3 ) . TURF ( 13) 

DIMENSION  MKASE( 100) .ZCOLAP( 100) .ZO3LAP(100) .ZFWV(IOO) .WFWV( 100) , 
1ZF03(100) .WF03( 100) .ZFCO2(100) ,WFCO2(100) ,ZTOTAL( 100) 

DIMENSION  RAT  101 (100) .RAT  102 (100) .RAT  103 ( 100) 

READ (50.60)  ( ( WATER ( I.J). J= 1.13). 1=1.21) 

READ (50.60) ( (C02( I . J ) . J=l » 13 > , I =1 .23 > 

READ (50 .60) ( ( OZONE ( I.J).J=l.l3).I=1.16) 

60  FORMAT( 13F6.3) 
READ(50.61) (U( I ).I=1.21 ) 
READ (50.61) (UUU( I ) .1=1,23) 
READ(50.61)   (UU( I ) .1=1.16) 

61  FORMAT ( 12F4.0.3F3.0.F2.0.7F3.0) 
READ(50.59)  ( RWATER ( K ) .K= 1 .21 ) 
READ (50 .59) ( RC02 ( K ) »K=1 . 2  3 ) 
READ (50.59) (ROZ(K) »K=1.16) 

59  FORMAT( 12F6.2) 

READ  (50.701)  ( COXMI S (  I  ) . I  =  1  .2 1  ) 
701  FORMAT(5F15.10) 

READ(50.703)  (UC (  I  )  .  I  =  1  .21  ) 
703  FORMAT( 16F4.2) 

DO  62  1=1.13 


42 


62  TT(  I  )=I*10-90 
WRITE(51.100) 

100  FORMAT! 1H1) 

WR1TE(51.63) ITT(I) .1=1.13) 

63  FORMAT  (4X.13F8.0) 
WRITE(51.100) 

WRITE  I  51.64)  (U( I ) »< WATER ( I. J), J=l, 13). 1=1. 21) 

64  F0RMAT(F5.1.13F8.3) 
WRITE<51»100) 

WRITE( 51.64) (UU( I ) .(OZONE! I. J), J= 1.13). 1=1. 16) 

WRITE(51.100> 

WRITE (51.64) (UUU(  I ). (C02(  I .J). J= 1.13) .1  =  1. 23) 

INTER=63 

KUT  =  0 

60NST=0.4838 
777  READ  (50.1000)  KZ 
1000  FORMAT ( 14) 

READ  (50.7000)  (COUNT 
7000  F0RMAT( 14) 

10  READ  (50.1)  (KASE<  I > .LEVEK I  I »P( I ) »T(I ) »Q( I ) .CQ( I ).I=1.KZ) 
1  FORMAT(2I3»F9.1.F5.0.F12.4,Fl0.3) 

WRITE  (51.100) 

WRITE  (51.2)  (KASE( I  ) .LEVEL ( I ) ,P ( I ) .Q< I ) ,CQ(  I  )  » I  =  1 » KZ  ) 
20FORMAT(1X.4HCASE.4X.5HLEVEL.4X,8HPRESSURE.6X.3HH20,7X.5HOZONE//(1X 
l.I4.I8.F14.4,F12,4,F11.3)) 

KT-  KZ  -  1 

3UMW=0. 

iUMO=0. 

SUMC=0. 

CONST  1  =  0. 4342 9448 19 

DO  2C  1=1. KJ 

PMEAN(  I  )  =  (P( I  )+P( 1+1) )/2. 

QMEAN(I)  =(Q( I )+Q(  1  +  1)  )/2. 

CQMEAN(I)  =  (CQ( I )+CQ( 1+1 ) )/2. 

PDIFF(  I )  =  Pi  I )-P( 1  +  1) 

SDIFFt I )  =  Q( I )-Q( 1+1) 
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CQDIFFI I )  =  CQ<  I )-CQ( 1  +  1) 

c     the  purpose  of  thes  sub-program  is  to  calculate  the  log  values  for 
c     the  reduced  optical  paths. 

i>UMW  =  SUMW+QMEAN(  I  )*PMEAN(  I  )*PDIFF(  I  )*<  1.00705E-6) 
SUMO=SUMO+CQMEAN( I )*PDIFF( I >*( l.E-5) 
bUMC=SUMC+PMEAN( I )*PDIFF( I ) * ( 2 • 4469E-4 ) 
SFW( I )=SUMW 
SFO( I )=SUMO 
SFC( I )=SUMC 
IF  (SFW(D)  56.57.58 
58  CFW( I )=CONSTl*LOGF(SFW( I ) ) 

56  CONTINUE 

IF  (SFO(I))  52.53.54 
54  CFO( I)«CONSTl*LOGF(SFO( I > ) 

52  CONTINUE 

IF  (SFC(I))  46.47.48 
48  CFC< I )«CONSTl*LOGF(SFC( I ) ) 

46  CONTINUE 
GO  TO  20 

57  CFW(I)*19. 
GO  TO  56 

53  CF0(I)=19. 
GO  TO  52 

47  CFC(I)=19. 
20  CONTINUE 

GWRITE151.4)  (I.T(I).P(I).Q(I) ,PMEAN( I ) ,PDIFF( I ) ,QMEAN( I ) .SFWU) . 
1SF0(  I  )  »SFC(  I  )  •  I  »1  »K.T  ) 
40FORMAT( 2H  I »4X . 1HT »6X . 1HP . 7X » 1HQ*4X. 5HPMEAN. 3X . 5HPDI FF , 3X . 
l5HQMEAN»6X»3HSFW.6X.3HSF0»6X.3HSFC//( I3.F6.0.F8.1»F7.3»2F8.1.4F9.4 
2)  ) 
WRITE(51.100) 

WRITE!  51  »9  J  (I.T(I  +  1)  ,CFW(  I)  »CFO(  I  )  ,CFC(  I  ).I  =  1»K,T) 
90FORMAT(3H   I.2X.6HT(  1  +  1 ) » 7X , 3HCFW. 17X.3HCFO. 17X.3HCFC/ / ( I  5. F 5.0. 
13E2  .6)  ) 

C      NOTE ANY  LOG  VALUE  DENOTED  AS  19  INDICATES  THAT  IT  HAS  BEEN 

C      FORCED  OUT  OF  RANGE  OF  THE   R   VALUE  TABLES  BECAUSE  OF  THE   U   VA- 
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L      LUE  bEINu  Zt.RO-   ALSO.  NOTE  THAT   U    VALUES  CAN  NOT  BE  NEGATIVE 
C      DUE  TO  THE  ARITHMETIC  AND  SEUUENTIAL  PROCESSES  OF  THE  PROGRAM. 

WRITE151.100) 

JUMP=1 

I^H  = 

M3  = 

K.C  = 
C      THE  PURPOSE  OF  THIS  SUB-PROGRAM  IS  TO  DETERMINE  THE  SPECIFIC 
C      HATER   R   VALUES  FOR  EACH  LEVEL  OF  EVERY  CASE- 
DO  13  l=K.OUNT»ia 

TU  =  T( 1+1)  -  273. 

wU=CFW( 1 ) 
ifa  DO  23  K.  =  l.i2 

IF(  TU-TT(K.)  )  23-40.50 
5C  IFITU-TT(K+1) )  40-43.23 
23  CONTINUE 

IF(JUMP)289.189.89 

39  WRITE<51.88)  KASE  (  I  )  ♦  T  U  +  l  )  »WU 

880FORMAT( 1X.8H  CASE  =  I4-4X.8H  TUH20  =F5.0.5X.11H  LOG  UH20  =F6.3.23 
1H  ARE  OUT  OF  TABLE  RANGE) 
133  IF(JUMP)107.26.13 

40  IT=K 

GO  TO  65 
43  IT  =  tC+l 

GO  TO  65 
lul  DO  39  L=1.20 

IF(WU-U(L> )39. 70.80 
80  IF  <WU-U(L+1))  70.73.39 
39  CONTINUE 

GO  TO  89 
70  IW  =  L 

74  IF( JUMP1120.110.90 
73  IW  =  L+l 

SO  TO  74 
900CALL  FFIN(WU.U( IW) »U(  IW+1) »TU.TT(IT).TT( IT  +  1).WATER(  IW.  IT)  . 
1 WATER ( IW.IT  +  1  )  .WATER  I  IW+1. IT) .WATER ( I W+l . I T  +  l ) . FUT ( I  )  ) 
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K.H  =  MH+1 

RH201KH)  =  FUT< I ) 

TEMPI (KH)=T( 1+1) 
13  CONTINUE 

WRITE(51»100) 

JUMP=JUMP-1 

DO  26  I=tCOUNT»KT 

TU=T<  I+D-273. 

WU«=CFO(  I) 

GO  TO  36 
189  WRITE(51»188)  KASE  (  I  )  »T  (  1  +  1 )  »WU 

18buFORMAT< 1X.8H  CASE  ■  I4»4X»8H  TU03  =  F5.0.5X.10H  LOG  U03  -F6.3.23H 
1ARE  OUT  OF  TABLE  RANGE) 

GO  TO  133 

102  DO  105  L=ltl5 
IFIWU-UU(L) )105»70.106 

106  1F(WU-UU(L+1) ) 70.73.105 
105  CONTINUE 

GO  TO  189 
110UCALL  FFOZ(WU.UU(  IW> .UU< IW  +  1 >  »TU.TT( IT) »TT< IT  +  1 ) »020NE( IW.IT  )  . 
lfc»ZONE( IW.IT  +  1  ) » OZONE ( IW+lt IT ) ,0ZONE( IW  +  1.  IT+1) *GUT( I ) ) 

K>0  =  K.0+1 

R03(K.0)=GUT(  I  ) 

TEMP2(K.O)  =  Tl  1  +  1) 
26  CONTINUE 

JUMP=JUMP-1 

DO  107  I=KOUNT.KT 

TU  =  T(  I  +  D-273. 

*U=CFC( I ) 

60  TO  36 
289  WRITE151.288)  KASE ( I ) t T ( 1+1 > tWU 

288  F0RMAT(1X,8H  CASE  =  I4.4X.8H  TUC02  =  F5.0t5X*llH  LOG  UC02  =F6.3.23 
1H  ARE  OUT  OF  TABLE  RANGE) 

GO  TO  133 

103  DO    108    L=l. 22 
IF(WU-UUU(L)  )     108.70.109 
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109  IF(WU-UUU(L+1 ) )70.73,108 
108  CONTINUt 

GO  TO  289 
12  00CALL  FFC02(WU.UUU( IW ) .UUU( IW  +  1) . TU .TT ( I T ) »TT ( IT+1),C02< IW.IT). 
1C02( IW. IT+1 ) »C02( IW+1 ♦ IT) »C02( IW+l.IT+1 ) .DUT( I ) ) 

K.C  =  K.C+1 

RRCO(XC)=DUT( I ) 

TEMP3(KC)  =  T(l+1) 
107  CONTINUE 

GO  TO  77 
6b  IF(JUMP)103.102.101 

c     this  sub-program  determines  the  flux  integral  for  water  vapor 

C      BETWEEN   -80  TO  -273  DEGREES  CENTRIGADE. 
77  WU  =  CFW(KLT) 
IND  =  3 
DO  150  1=1.20 
IF  (WU  -  U( I ) )  150.170.180 

180  IF(WU-U( 1+1 ) )  170.173.150 
150  CONTINUE 

L70  IW=I 

174  GO  TO  (220.210. 190)IND 

173  IW=I+1 

GO  TO  174 

190  CALL  RRWAT(WU»U( IW) »U( IW+1) »RWATER( IW) .RWATER( IW+1) »FF) 
WRITE  (51.191)  FF 

191  F0RMAT(1X»6H  FF  =  E17.9) 

C     THIS  SUBPROGRAM  DETERMINES  THE  FLUX  INTEGRAL  FOR  OZONE  BETWEEN 
C      -80  TO  -273  DEGREES  CENTRIGADE. 

IND=IND-1 

WU=CFO(KT ) 

DO  160  1=1.15 

IF(WU-UU( I ) ) 160.170. 181 

181  1F(WU-UU( 1+1) ) 170.173.160 
160  CONTINUt 

210  CALL  RROZ(WU.UU( IW).UU( IW+1) .R02( IW) »ROZ( IW+1) .GG) 
WRITE  (51.211)  GG 
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211  F0RMAT(1X.6H  GG  =  E17.9) 
C      THIS  SUB-PROGRAM  DETERMINES  THE  FLUX  INTEGRAL  FOR  C02  between 
C      -8U  TO  -273  DEGREES  CENTRIGADE. 
IND=IND-1 
WU  =  CFC(K,T  ) 
DO  200  1=1.22 
IF(WU-UUU( I ) )  200.170.182 
182  IF(WU-UUU( 1+1 ) )  170.173.200 
200  CONTINUE 

220  CALL  RRC02(WU.UUU( IW) »UUU( IW+1) »RC02( IW) »RC02( IW+1) .DD) 
WRITE  (51.221)  DD 

221  FORMAT! IX. 6H  DD  =  E17.9) 

C      THIS  SUB-PROGRAM  INTEGRATES  THE  INTEGRAL   RDT   BETWEEN  T   INITIAL 
C      TO   T   FINAL  IN  DEGREES  CENTRIGADE. 

IND=IND-1 

AREA1=0. 

AREA2=0. 

AREA3=0. 

TOTAL=0. 

fcHH=KH-l 

DO  400  I  =  1»K.HH 

Y»TEMP1(  I  )  -  TEMPK  1  +  1  ) 
40  0  AREA1«AREA1+(RH20( I )+RH20( 1+1 ) ) /2.*Y 

WRITE  (51.413)  AREA1 

413  FORMAT!  IX. 10H  AREA1   =  E17.9) 
<00=K0-1 

DO  401  I^l.KOO 
Y=TEMP2( I )  -  TEMP2( 1  +  1 ) 
401  AREA2  =  AREA2+(R03(  I  )+R03(  1+1) >/2«*Y 
WRITE  (51.414)  AREA2 

414  FORMAT( 1X.10H  AREA2   =  E17.9) 
K,CC  =  KC-1 

DO  402  1=1. KCC 
Y*TEMP3( I )  -  TEMP3( 1  +  1  ) 
40  2  AREA3«=AREA3+(RRCO(  I  )+RRCO(  1  +  1)  )/2.*Y 
WRITE  (51.415)  AREA3 
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h\b    FORMAK  1X.10H    AREA3       =    E17.9) 

c    this  sub-program  integrates  the  integral  rdt  between  t  final 

C      TO   80  DEGREES  CENTRIGADE. 

MC  = 

MO= 

MW  = 

KUMP*1 

WU-CFWIKLT) 

DO  900  L=1.13 

DO  903  J=1.20 

IF<WU-U(J>)  903,906.909 
909  IF(WU-U< J+l ) )906»973.903 
903  CONTINUE 

GO  TO  900 
906  1W=J 

974  IF(KUMP)  983.982.981 
973  IW=J+1 

GO  TO  974 

981  CALL  FFFNIWU.UI  1W )  ,U(  IW  +  l  ) > WATER (  IW.L) .WATER (  IW+l.L)  .ROPE(L)  ) 
MW=MW+1 

TTW(MW)  *  TT(L) 

RW(MW)  =  ROPE(L) 
900  CONTINUE 

K.UMP=KUMP-1 

WU=CFO(KT> 

DO  915  L=1.13 

DO  918   J«1.15 

IF(WU-UU(J))  918.906.924 
924  IF(WU-UU< J+l ) )  906.973.918 
918  CONTINUE 

GO  TO  915 

982  CALL  FFFN(WU.UU( IW).UU< IW+l) .OZONE ( IW.L) .OZONE < I W+l . L ) .SOPE ( L ) ) 
MO=MO+l 

TTO(MO)=TT(L) 
RO(MO)  =  SOPE(L) 
915  CONTINUE 
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K>UMP  =  KUMP-1 

wu=CFC<ia) 

DO  930  L-1.13 

DO  933  J»l»22 

IF <*U-UUU< J) )  933.906.936 
936  IF(WU-UUU( J+l ) )  906.973.933 
933  CONTINUE 

GO  TO  930 
98  3  CALL  FFFNCWU.UUUJ  IW) .UUU( IW+1) ,C02 ( IW.L) »C02( IW+l.L) .TOPE(L) ) 

MC=MC+1 

TTC(MC)  =  TT(L) 

RC(MC)  =  TOPE(L) 
930  CONTINUE 

MUMP=1 

TU  =  TEMPKK.H)    -273. 

KL  =  KH 

i^M  =  KO 

ttN  =  KC 
830    DO    800    K=1.12 

IFITU-TT(K))    800.803.806 
806    IF(TU-TT(K+1) )     803.873.800 
800    CONTINUE 

IF(MUMP)    815.812.809 
809    *RITE(51»818)     TU 
818    FORMATdOH    TF    H20    =    F5.0.3X.22H    IS    OUT    OF    TABLE    RANGE//) 

M_  =  K.L-1 

TU  =  TEMPKKL)-273. 

GO    TO    830 
803     IT  =  K. 

874  IF(MUMP)  883.882.881 
873  IT=fc+l 

GO  TO  874 
881  CALL  F0RGE1(TU.TTW(  IT)  »RW(  IT)  ,RH20(ICL)  .AREA11) 

ADD1=0. 

IM= I T-l 

DO  824  L=1.IM 
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CALL  F0RGE2(RW(L) .RW(L-H) » TTW( L ) »TTW < L+ 1 ) ,RURF ( L ) ) 
824  CONTINUE 

DO  827  L=1»IM 

827  add1=add1+rurf(l) 

this  sub-program  determines  the  total  flux  for  water  vapor  with- 

0ut  c02  and  ozone  overlap  corrections. 
845  area1  =  area1  +  add1  +  area11 

write  (51*410)  add1,area11,area1 

410  f0rmat(1x»8h  add1  *  e17.9»5x»10h  area11  =  e17.9 ♦ 5x» 14ht0tal  area1 

1=    E17.9) 

MUMP=MUMP-1 

TU=TEMP2(KO)    -273. 

GO    TO    8  30 
812    *RITE(51.833)     TU 
833    FORMATUOH    TF       03    =    F5.0.3X.22H    IS    OUT    OF    TABLE    RANGE//) 

KM=KM-1 

TU=TEMP2(KM)    -273. 

GO  TO  830 
882  CALL  FORGEUTU.TTO<  IT)  .RO(  IT)  , R03 < KM ) ♦ AREA22 ) 

ADD2=0. 

IM=IT-1 

DO  836  L=1.IM 

CALL  FORGE2(RO(L) .ROJL+1) • TTO( L ) *TTO( L+l ) »SURF ( L )  ) 
836  CONTINUE 

BO  839  L=1»IM 
839  ADD2=ADD2+SURF(L) 

THIS  SUB-PROGRAM  DETERMINES  THE  TOTAL  FLUX  FOR  OZONE. 
848  AREA2«AREA2  +  ADD2  +AREA22 

WRITE  (51*411)  ADD2.AREA22»AREA2 

411  FORMAT! IX. 8H  ADD2  =  E17.9.5X»10H  AREA22  =  El  7.9. 5X» 14HT0TAL  AREA2 
1=  E17.9) 

MUMP=MUMP-1 

TU=TEMP3(KC)  -273. 

GO  TO  830 
815  WRITE(51.842)  TU 
842  FORMAT (10H  TF  C02  =  F5.0.3X.22H  IS  OUT  OF  TABLE  RANGE) 


51 


KN=KN-1 

TU*TEMP3<KN>  -273. 

60  TO  830 
883  CALL  F0RGE1<TU»TTC(IT) tRC( IT) ,RRC0<KN) tAREA33) 

ADD3«0. 

IM=IT-1 

DO  851  L=1»IM 

CALL  F0RGE2(RC(L) »RC(L+1) » TTC ( L ) »TTC ( L+l ) .TURF(L) ) 
851  CONTINUE 

DO  854  L=1»IM 

854  add3=add3+turf(l) 

this  sub-program  determines  the  total  flux  for  carbon  dioxide. 
857  area3=area3+add3+area33 

write  (51*412)  add3»area33.area3 
412  f0rmat(1x,8h  add3  =  e17.9»5x»10h  area33  =  e17.9 , 5x , 14ht0tal  area3 

1=  E17.9) 

THIS  SUB-PROGRAM  DETERMINES  THE  C02  and  ozone  overlap  CORRECTIONS. 

JIP  =  3 

DO  450  1=1,20 

WU-CFW(KT ) 

IF(WU-UC(I)>  450»470,480 
480  IF(WU-UC( 1+1) )  470»473»450 
45u  CONTINUE 

MRITE(51*451)WU 

451  F0RMAT(17H  LOG  UF  TAUC02  =  F9.4»3X*22H  IS  OUT  OF  TABLE  RANGE//) 
RCOLAP«0. 

GO  TO  452 
470  IW=I 
474  GO  TO  (490t491»492 )JIP 

490  CONTINUE 

491  CONTINUE 
473  IW=I  +  1 

60  TO  474 
4920CALL  AOVLAP I WU»UC ( I W ) »UC< Ito+1) ,COXMIS( IW) .COXMISi IW+1 ) ♦COTAUF(KT) » 
1AUD3»AREA33.DD.RC0LAP) 

452  JIP=JIP-1 
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00  886  I=K0UNT»KT 
886  ©3TAUFU)  =  EXPF(-0.1167  *  SFW(D) 
R03LAP=03TAUF ( KT ) * ( ADD2+AREA22+GG) 

MP  =  3 

DO  549  J-KOUNT.KT 

WU=CFW( J) 

DO  550  1=1.20 

IF(WU-UC(IU  550.570,580 
580  IF <WU-UC( 1+1 ) )  570.573*550 
550  CONTINUE 

WRITE(51.453)  ICASE(J).WU 
4530FORMAT( 1X.7H  CASE  =I4.3X.16H  LOG  UC02  TAU  =»  F9.4.22H  IS  OUT  OF  TAB 
ItE  RANGE) 

60  TO  549 
570  IW  =  I 
574   GO  TO  (590,591.592)KIP 

590  CONTINUE 
573  IW=I+1 

GO  TO  574 
5920CALL  BOVLAP ( WU.UC< IW) »UC( IW+1) .COXMISl IW)»COXMISl IW+1). 
1C0TAUF(  J)  .DUT(J).TU)  »T(J  +  1)  .C02LAP ( J ) .DUT ( J+l )  ) 

KV=KV+1 

C02LAP(KV)  =  C02LAP(J) 
549  CONTINUE 

KIP  =  KIP  -  1 

TCOLAP  =  RCOLAP 

DO  449  J=1.KV 
449  TCOLAP  *  TCOLAP  +  C02LAP(J) 

DO  454  J=KOUNT.KT 

591  CALL  BBVLAP(03TAUF(J) »GUT(J) »GUT( J  1 ) .T ( J ) .T ( J+l ) .03RLAP ( J  )  ) 
454  CONTINUE 

T03LAP  «=  R03LAP 
DO  569  J*(COUNT.KT 
569  T03LAP  =  T03LAP  +  03RLAP(J) 

THIS  SUB-PROGRAM  DETERMINES  THE  TOTAL  DOWNWARD  FLUX  FOR  EVERY 


53 


SOUNDING. 

CALL  ASUMF ( ARE Al . ARE A2 .AREA3 »FF .GG.DD. TCOLAP .T03LAP . TOTAL) 

WRI TEC  51*79) ( AREA 1 .AREA2 .AREA3 .FF »GG*DD» TCOLAP » T03LAP . TOTAL ) 
790FORMAT(2X»5HAREA1»3X,5HAREA2»3X»5HAREA3,6X»2HFF,7X,2HGG»6X»2HDD» 
15X.6HTCOLAP»2X,6HTO3LAP»2X.10hTOTAL  FLUX// ( IX, 7E 17.9 ) ) 

FWV  =  AREAlfFF-TCOLAP/2.-T03LAP/2. 

F03  =  AREA2+GG-T03LAP/2. 

FC02  =  AREA3+DD-TC0LAP/2. 

THE  FOLLOWING  VALUES  HAVE  UNITS   CAL/CM  SQUARE/DAY 

WRITE  (51.201 )  AREA1 fFF.AREA2,GG,AREA3» DO. TCOLAP, T03 LAP. FWV.F03* 
1FC02. TOTAL 
2010FORMAT(1X,10H  AREA1   =  E17.9//1X,1  H  FF      =  E 17.9//1X, 10H  AREA2 

1  *  E17.9//1X,10H  GG      =  E17 . 9//1X , 10H  AREA3   =  E17.9// IX  ,  10H  DD 

2  =  E17.9//1X,10H  TCOLAP  =  E17.9// IX , 10H  T03LAP  =  El  7.9// IX, 10H 
3FWV     =  E17.9//1X,10H  F03     »  E17.9// IX ,  10H  FC02    *  E17.9///1X, 
423H  TOTAL  DOWNWARD  FLUX  =  E17.9) 

XAREA1  =  AREA1*C0NST 
XFF  =  FF*CONST 
XAREA2  =  AREA2*C0NST 
XGG  =  GG*CONST 
XAREA3  «=  AREA3*C0NST 
XDD  =  DD*CONST 
XCOLAP  =  TCOLAP*CONST 
X03LAP  =  T03LAP*CONST 
XFWV  =  FWV*CONST 
XF03  «=  F03*CONST 
XFC02  =  FC02*CONST 
XTOTAL  =  TOTAL*CONST 

THE  following  values  have  units  watts/meters  square 

WRITE(51.1011)  K.ASE(ICT) 
1011  FORMATC 1X.7HCASE  =  14) 

WRITE! 51*10  01 )  XAR EA 1.  XFF »XAREA2*XGG»XAREA3» XDD. XCOLAP. X03LAP. XFWV 
1  .XF03.XFC02.XT0TAL 
lCul  FORMATC 1X.10H  NAREA1  =  E17.9//1X.1  H  NFF     =  E17.9//1X. 1  OH  NAREA2 

1  =  E17.9//1X.10H  NGG     -    El  7. 9// IX , 10H  NAREA3  =  El  7.9// IX , 10H  NDD 

2  =  E17.9//1X,10H  NCOLAP  =  E  1  7. 9// IX ♦ 10H  N03LAP  =  E17. 9/ / IX, 10H 
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3NFWV    =  E17.9//lXtlOH  NF03    »  El  7. 9// IX . 10H  NFC02   «  E17.9//lXt 
<t23H  TOTAL  DOWNWARD  FLUX  =  E17.9///) 

INTER=INTER-1 

fvUT  =  K.UT  +  1 

MKASE(ICUT)     =    KASEOCT) 

ZCOLAP(KUT)  =  XCOLAP 

Z03LAP(KUT)  =    X03LAP 

ZFWV(KUT)  =  XFWV 

ZF03(KUT)  =  XF03 

ZFC02(KUT)     =    XFC02 

ZTOTAL(KUT)     =    XTOTAL 

WFWVUUT)    =    XAREA1    +    XFF 

WF03(KUT)     =    XAREA2    +    XGG 

WFC02(K.UT)     =    XAREA3    +    XDD 

IF    (INTER)    888,999.777 
999    WRITE    (51»100) 

YKUT    =    (CUT 

ADD1    =    0.0 

ADD2    =    0.0 

ADD3    =    0.0 


AVG11  = 

0.0 

AVG22  = 

0.0 

AVG33  = 

0.0 

AVG44  * 

0.0 

AVG55  = 

0.0 

AVG66  = 

0.0 

AVG77  = 

0.0 

AVG88  = 

0.0 

AVG99  = 

0.0 

DO  8050 

I=1»KUT 

AVG11  = 

AVG11 

+ 

ZC0LAP1 I ) 

AVG22  = 

AVG22 

+ 

Z03LAP( I) 

AVG33  = 

AVG33 

+ 

ZFWV( I ) 

AVG4-4  = 

AVG44 

+ 

WFWV( I ) 

AVG55  = 

AVG55 

+ 

ZF03( I ) 

AVG66  = 

AVG66 

+ 

WF03( I ) 
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AVG77  *  AVG77  +  ZFC02 <  I  ) 
AVG88  =  AVG88  +  WFC02 (  I  ) 
AVG99  =  AVG99  +  ZTOTAH  I) 
8050  CONTINUE 

AVG1  «  AVG11/YXUT 

AVG2  =  AVG22/YKUT 

AVG3  =  AVG33/YK.UT 

AVG4  =  AVG44/YKUT 

AVG5  =  AVG55/YKUT 

AVG6  =  AVG66/YKUT 

AVG7  =  AVG77/YKUT 

AVG8  =  AVG88/YKUT 

AVG9  =  AVG99/YK.UT 

DO  1111  I»1»KUT 

RATIOKI)  =  ZFWV( I )/ZTOTAL< I ) 

RAT  102(1)  =  ZF03(  I )/ZTOTAL( I ) 

RATI03U)  =  ZFC02I I )/ZTOTAL( I ) 

1111  CONTINUE 

BO  1112  1  =  1. KUT 
A001  =  AD01  +  RATIOK  I  ) 
A0D2  =  ADD2  +  RATI  02 (  I  ) 
ADD3  =  ADD3  +  RATI031  I  ) 

1112  CONTINUE 

A0D11  »  <ADD1/YKUT)*100. 

ADD22  =  ( ADD2/YKUT )*100. 

A0033  =  (ADD3/YKUT )*100. 

WRITE(51»1022>  (MKASE( I )*ZCOLAP< I ) »Z03LAP( I ) *WFWV<  I  )  *ZFWV<  I  )  , 
lwF03(  I  )  tZF03(  I  )  »WFC02(  I  )*ZFC02(  I)  t ZTOTAH  I)  *I=1»KUT) 
1022  FGRMAT< 1X.85HCASE   H20-C02   H20-03   UNCOR-    CORRECT-   UNCOR-    CO 
1RRECT-   UNCOR-    CORRECT-   T0TAL/1X.88H       OVERLAP   OVERLAP   REC 
2TED    ED  WV     RECTED     ED  03     RECTED     ED  C02    D0WNWARD/1X.8 
3AH  WV  FLUX    FLUX      03  FLUX    FLUX      C02 

A  FuUX   FLUX      FLUX/ ( 1X» I4*2F8. 3* F10. 3 » 3F9. 3 . 3F10.3  )  ) 

KRITE(51*1031)  AVGl*AVG2.AVG3.AVG<f»AVG5»AVG6*AVG7»AVG8.AVG9 
1031  FORMAT ( IX . ^HMEAN/IX* 5H VALUE ♦ F7.3 »F8 . 3 *F10. 3 » 3F9. 3 » 3F10.  3 ) 

WR1TE( 51*1032)  ADD11 .ADD22 »ADD33 


56 


1U32  FURMAK 1X»31HAVG  VALUE  CONTRIBUTION  ON  TOTAL/1X.32H   DOWNWARD  FLUX 

1  IN  PERCENT .F6.2.F16.2.F20.2 ) 

B88  END 
C      THIS  SUbROUTINE  EXECUTES  FORWARD  AND  BACKWARD  INTERPOLATION  TO  OBTAIN  THE 
C      SPECIFIC  ELASSER  WATER  R   VALUE  FOR  EVERY  REDUCED  OPTICAL  PATH/TEMPERA- 
C      TURE  PAIR. 

SUBROUTINE  FFI N ( U.U1 »U2 »T . Tl »T2 »FU1T1 .FU1T2 .FU2T1 .FU2T2 .F ) 

FFUT1=FU1T1+(FU1T2-FU1T1)*ABSF( ( T-Tl ) / ( T2-T1 ) ) 

FFUT2*FU2T1+(FU2T2-FU1T1)*ABSF(  ( T-Tl  )  / ( T2-T1  )  ) 

F1=FFUT1+(FFUT2-FFUT1)*ABSF( IU-U1 ) / (U2-U1 ) ) 

FFUT3  =  FU1T2-(FU1T2-FU1T1)#ABSF((T2-T)/(T2-T1) ) 

FFUT^=FU2T2-(FU2T2-FU2T1)*ABSF( < T2  T)/(T2-T1 ) ) 

F2=    FFUT3+(FFUT4-FFUT3)*ABSF( ( U-Ul ) / (U2-U1 ) ) 

F=(F1    +    F2) /2. 

END 
USUBROUTINE  FFOZ (U.U1 »U2 »T »T1 »T2 »GU1T1 »GU1T2 »GU2T 1 »GU2T2 »G) 

GGUT1  =  GU1T1+(GU1T2-GU1T1)*ABSF( ( T-T 1 ) / ( T2-T1  )  ) 

GGUT2  =  GU2T1+(GU2T2-GU2T1 )*ABSF( ( T-Tl ) / ( T2-T1  )  ) 

G1=GGUT1+(GGUT2-GGUT1)*ABSF( (U-Ul ) / ( U2-U1 ) ) 

GGUT3=GU1T2-(GU1T2-GU1T1)*ABSF( <T2  T)/(T2-T1 ) ) 

GGUT4=GU2T2-(GU2T2-GU2T1)«ABSF( (T2  T)/(T2-T1 ) ) 

G2=GGUT3+(GGUT4-GGUT3)*ABSF( (U-Ul ) / (U2-U1 ) ) 

G=(Gl+G2)/2. 

END 
05UBROUTINE    FFC02 (U »U1 »U2 . T .Tl ,T2 »DU1T1 .DU1T2 .DU2T 1 ,DU2T2 »D> 

DDUT1=DU1T1+(DU1T2-DU1T1)*ABSF( ( T-Tl ) / ( T2-T1 ) ) 

DDUT2    =    DU2T1    +    (DU2T2    -    DU2T1)    *    ABSFUT-T1)    /    (T2    -    Tl  )  ) 

D1=DDUT1+(DDUT2-DDUT1)*ABSF( (U-Ul ) / (U2-U1 ) ) 

DDUT3=DU1T2-(DU1T2-DU1T1)*ABSF( (T2    T)/(T2-T1) ) 

DDUT4=    DU2T2-(DU2T2-DU2T1)*ABSF( ( T2-T ) / ( T2-T 1 ) ) 

D2    =DDUT3+(DDUT4-DDUT3)*ABSF( ( U-Ul ) / (U2-U1 ) ) 

D=(Dl+D2)/2. 

END 

SUBROUTINE  RROZ ( U.U1 »U2 .ROU1 .ROU2 »GG) 

GG=ROUl+(ROU2-ROUl )*ABSF( (U-Ul ) /(U2-U1 > ) 

END 
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SUBROUTINE  RRWAT ( U.U 1 .U2 .RWU1 .RWU2 .FF ) 
FF=RWU1+(RWU2-RWU1)*ABSF( ( U-Ul ) / ( U2-U1 > ) 

END 

subroutine  rrco2<u.ui.u2.rcui,rcu2.dd> 

DD=RCU1+(RCU2-RCU1 )*ABSF( (U-Ul )/(U2-Ul ) ) 

END 

SUBROUTINE  ASUMF ( A»B »C t D. E »F tX • Y .Z ) 

Z  =  A+B+C+D+E+F-X-Y 

END 

SUBROUTINE  B0VLAP(U.U1»U2.R1»R2»ZZ.D1.T1»T2.XX.D2> 

11=    Rl  +  (R2-R1)*ABSF( (U-U1)/(U2-U1) ) 

XX=ZZ*(Dl+D2)/2.*( T1-T2) 

END 

SUBROUTINE  AOVLAP ( U.U1 tU2 .Rl »R2 »Z . A.B.DD .X ) 

Z=R1+(R2-R1)*ABSF(  (U-Ul  )/(U2~Ul>  ) 

X=Z*(A+B+DD) 

END 

SUBROUTINE  BBVLAP(A.R1.R2.T1»T2»Z) 

Z=  A*  ((Rl  +  R2)/2.  *  (T1-T2)) 

END 

SUBROUTINE  FFFN ( U.U 1 .U2 »R1  .R2 , XX ) 

XX=R1+(R2-R1)*ABSF( ( U-Ul ) / ( U2-U1 ) ) 

END 

SUBROUTINE  FORGE1 ( TF , T 1 »R1 ,RF ,X ) 

X=(R1+RF)/2.*(TF-TD 

END 

SUBROUTINE  F0RGE2(R1»R2.T1.T2,Y) 

Y  =  ABSF(  (R1  +  R2)/2.*(T2-T1)  ) 

END 

END 

FINIS 
■EXECUTE. 


606 


58 


Appendix  I  B.   Modification  Instructions  to  Program  RADFLUX  (FORTRAN  63) 

1.  Each  sounding  is  subdivided  into  levels  starting  from  the  interface 
at  level  I  =  1  to  the  top  of  the  atmosphere  (0.1  mb)  at  level  I  =  KZ. 

KZ  is  the  number  of  the  last  level  in  the  sounding. 

i 

2.  The  following  must  be  completed  in  order  to  determine  F,  from  any 

preselected  reference  level  different  from  the  interface  level: 

a.   Note  the  count  or  number  of  the  pre-selected  level  in  the  sound- 
ing and  introduce  this  number  as  shown  below  in  the  example . 
1.   Suppose  selected  reference  level  is  7  out  of  30  levels  in 
a  sounding,  hence,  punch  this  number  on  a  data  card  with 
the  unit  digit  in  column  4.   This  reference  level  card  is 
the  second  data  card  in  each  of  the  case  atmospheres.   Do 
this  procedure  anytime  the  reference  level  changes  for  any 
case-atmosphere . 

3.  If  a  level  is  added  or  subtracted  in  a  sounding  listing,  the  number 
on  the  first  data  card  of  the  respective  case  atmosphere  must  be  changed 
accordingly.   For  example: 

a.  Suppose  that  there  are  29  levels  in  a  sounding,  thus,  the  first 
data  card  in  the  sounding  must  read  29;  however,  suppose  that  1 
level  is  added  to  the  sounding,  thence,  the  first  data  card  must 
be  changed  to  read  30. 

b.  One  must  place  the  unit  digit  of  the  specific  number  in  column 
4  of  the  data  card. 

4.  Finally,  the  INTER  (constant)  card  must  replaced  when  the  number 
of  case  atmospheres  in  the  study  changes.   For  example: 

a.   Presently  there  are  63  case  atmospheres  as  data  input  to  program 
RADFLUX,  thus,  the  INTER  card  reads  "INTER  =  63." 
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b.   Suppose  the  data  input  changes  to  100  case  atmospheres;  thence, 
the  INTER  card  must  be  changed  to  read  "INTER  =  100." 
5.   Program  RADFLUX  was  written  for  FORTRAN  63,  however,  this  program 
can  be  made  operational  in  any  FORTRAN  computer  language  provided 
that  minor  modifications  to  the  control  cards  are  made.   One  must 
check  with  the  respective  computer  language  manuals  for  further 
information. 
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Appendix  I,   Section  C.      Computer  Program  SUMFLUX   (FORTRAN  63) 


-C00P»B0X    T.    TUPAZ    J *S/ 1S/2S »20 . 20000 . 
-FTN.L*E. 

PROGRAM    SUMFLUX 

c    the  purpose  of  this  program  is  to  integrate  radiation  intensity 

C      CASE. 

C      BETWEEN  0  AND  90  DEGREES  TO  OBTAIN  TOTAL  UPWARD  FLUX  FOR  EVERY 

ODIMENSION  KASE<7) »PRES(7) »THETA(7) »X(7) »X2(7) ,UH2(7)  »U03(7)*X4(  7) ♦ 
1TEMP(7) .TOTINT(7) »CH4INT<7) 

DIMENSION  KOUNT(IOO) .PRESS (100) » TEMPER ( 100),CHFLUX( 100)  • 
1BBFLUX( 100) 

INTER  =  62 

K.UT  =  0 

CONST=0.56687E-7 

PI  =  3.14159 

RADCON  =  0.017453 

KZ  -  6 

THETA(6)=    90.0 
7  7  70READ    <50.2)     K.ASE<1  )»PRES(1  )»ThETA(  1)  »X(  1)  »X2 ( 1 ) »UH2 ( 1 ) »U03 < 1 ) »X4( 1 
1 » .TEMP( 1) 

2  FORMAT ( I3.F9.1.F8.0»F10.3.2F7.3»F6.3.F10.3»F7.1) 

•a  =  Ki-i 

OREAD  (50»3)  ( KASE ( I ) »THETA ( I ) , X ( I ) »X2 < I ) »X4< I ) ♦ I =2 »KT ) 

3  FORMAT ( I3»9X.F8.0»F10.3»F7.3tl3X.F10.3) 
SUM1=0. 

SUM2=0. 

DO  100  I»ltKT 

CALL  FORD1 (XI 3) »X ( 4 ) »X ( 5 ) »X4 ( 3 )  »X4 ( 4 ) »X4 ( 5 ) »X ( 6 ) »X4 ( 6  )  ) 

CALL  FORD2(X( I ) »X ( I  +  1 ) »X4 ( I ) »XA ( 1  +  1 ) »THETA( I ) ,THETA( 1  +  1 ) »TOT INT (  I  ) 
1«CH4INT(I) tRADCON) 

WRITE (51. 5)  KASE( I ). I . TOT  I  NT ( I ) »CH4INT(  I ) 
50FORMAT( 1X.8H  CASE  =  I5.5X.5H  I  =  I4»5XtlOH  TOTINT  =  E17.9.5X.10H  C 
1H4INT  *  E17.9//) 

SUM1=SUM1+T0TINT( I ) 

SUM2=SUM2+CH4INT( I ) 
100  CONTINUE 

WRITE(51»4)  <KASE(  1)  tPRES(l)  »TEMP(  1  )  ♦(  UTHETAi  I  )  »  X  (  I  )  »X4  (  I  )  .  I  =  1  ,KZ 
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1 » ) 

40FORMAT( 1X.8H  CASE  =  I5.5X.12H  PRESSURE  ■  F5.0.5X.8H  TEMP  =  F5.0// 
K1X.5H  I  =  I4.5X.10H  THETA1  *  F5.1.5X.19H  TOTAL  INTENSITY  =  F6.3. 
25X.18H  CH  4  INTENSITY  =  F6.3//)) 

FLUX1«SUM1*PI 

FLUX2*SUM2*PI 

WRITE(51.6)  KASEd  )  .FLUX1.FLUX2 
6  F0RMAT(1X.8H  CASE  ■  I5.5X.9H  FLUX1  =  E17.9.5X.9H  FLUX2  =  E17.9///) 

CALL  FORD3(CONST»TEMP( 1 ) .FLUX3) 

WRITE! 51.300)  FLUX3 
300  F0RMAT(1X.9H  FLUX3  ■  E17.9///) 

PRESLN  =  LOGF<PRES( 1 ) ) 

INTER  =  INTER-1 

KUT  =  (CUT  +  1 

KOUNT(KUT)  =  KASE<1  ) 

PRESS(KUT)  *  PRES( 1) 

TEMPER! KUT)  ■  TEMP(l) 

CHFLUX(KUT)  =  FLUX2 

BBFLUX(KUT)  =  FLUX3 

IF  (INTER)  888*999.777 
999  WRITE(51.200) 
200  FORMATt 1H1) 

WRITE (51.400)  <K0UNT(  I ) .PRESS ( I ) .TEMPER (  I ),CHFLUX( I ) .BBFLUX(  I).  1  =  1 
l.KUT) 
400  FORMAT( 1X.53HCASE    INTERFACE    INTERFACE     CHANNEL  4        BLACK 
1/1X.52H        PRESSURE    TEMPERATURE      FLUX  B0DY/49X. 

24HFLUX//( 1X.I4.F8.0.F13.0.F15.3.F14.3)) 
888  END 

SUBROUTINE  FORD1 ( A »B .C .D. E »F »G .H ) 

G  «=  C  +  0.6866*A  -1.8649*B  +1.1783*C 

H  =  F+0.6866»D  -1.8649*E  +1.1783*F 

END 

SUBROUTINE  F0RD2 ( A .B »C »D. E »F »X ♦ Y»Z ) 

X=(A+B)/2.*( (SINF(F*Z  )  )»*2  -  ( S I NF ( E*Z ) ) **2 ) 

Y  =  (C+D)/2.*( (SINF(F*Z) )**2  -  ( SI NF ( E*Z ) ) **2  ) 

END 
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SUBROUTINE  FORD3(A.B,C) 

OA*B**4 

END 

END 

FINIS 
■EXECUTE. 
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